First molecular detection and genetic diversity of Hepatozoon sp. (Apicomplexa) and Brugia sp. (Nematoda) in a crocodile monitor in Nakhon Pathom, Thailand

The crocodile monitor (Varanus salvator) is the most common monitor lizard in Thailand. Based on habitat and food, they have the potential to transmit zoonoses, with a high possibility of infecting ectoparasites and endoparasites. Diseases that could infect crocodile monitors and be transmitted to other animals, including humans. This research aims to identify and evaluate the phylogenetic relationships of Hepatozoon sp. and sheathed microfilaria in crocodile monitors. The phylogenetic analyses of Hepatozoon, based on 18S rRNA, and sheathed microfilaria, based on the COX1 gene, revealed that the Hepatozoon sp. were grouped with H. caimani, while sheathed microfilaria were grouped together with B. timori. This study provides insights into the genetic diversity and host-parasite interactions of hemoparasites in crocodile monitors in Thailand.

www.nature.com/scientificreports/ the morphotype, and thus, it was classified as an undescribed species.The examination of blood smears revealed the presence of mature gamonts, as depicted in Figs.1A,B.
In the blood smear, only mature gamonts were identified.These mature gamonts were found within a parasitophorous vacuole (PV) and had a rounded shape at both ends.The cytoplasm of the gamonts displayed a bluish-purple stain.The elongated nucleus exhibited purple-stained chromatin and occupied nearly half of the surface area of the parasite.Gamonts measured (mean ± standard deviation) 11.32 ± 0.886 in length and 4.25 ± 0.621 in width (n = 30) (Table 1).

Molecular amplification, sequencing, and similarity of Hepatozoon 18S rRNA and Brugia COX1 sequences
One captured crocodile monitor was observed with both hemogregarine and filarial infections (Fig. 1A,B).The polymerase chain reaction (PCR) products were sequenced for the 18S rRNA gene (hemogregarine) and COX1     3.

Phylogenetic analysis of Hepatozoon 18S rRNA and Brugia COX1 sequences
The phylogenetic tree of Hepatozoon sp.18S rRNA gene sequences demonstrated clustering in a monophyletic group together with sequences of H. caimani (KU495923) that were recently detected in caiman crocodiles in Brazil (Fig. 2).The phylogenetic tree was comprised of eight branches.and O. japonica (Fig. 3).The amplified filarial sequence was clustered with sequences of B. timori (AP017686 and KP760173), recently detected in humans in Japan and Indonesia.Nematoda sp. and Setaria sp. were used as the out-groups (Fig. 3).Moreover, the reliability of bootstrap frequencies and Bayesian posterior probabilities of all phylogenies are displayed with the highest values on each branch.

Haplotype diversity
Nucleotide polymorphisms and DNA divergence between the sequences obtained in the present study and Gen-Bank sequences were analyzed.Nucleotide polymorphism analysis of Hepatozoon sp.18S rRNA and Brugia sp.COX1 sequences revealed 28 and 24 haplotypes, respectively (Table 4).The haplotype networks of these genes were obtained from the Templeton, Crandall, and Sing (TCS) Network tool (Figs. 4, 5).For the Hepatozoon sp.18S rRNA gene, of the 28 haplotypes, haplotype 1 was detected in the crocodile monitor from the Nakhon Pathom provinces and in H. caimani in the caiman crocodile from Brazil.Haplotypes 1-23 were found in reptiles, rodents, and amphibians, while the remaining haplotypes were found in tick, canine, feline, and avian hosts from a range of countries (Table 3, Fig. 4).The haplotype network of Brugia sp.COX1 gene demonstrated that the sequence from the Thailand crocodile monitor was detected in haplotype 1, while the B. timori sequences from Japan and Indonesia were detected in haplotype 2 (Table 3, Fig. 5).

Discussion
Until now, the detection of Hepatozoon infections in crocodile monitors has been only reported in Bangkok, Thailand, by means of microscopy screening of blood samples 17 .Morphologic and morphometric studies could enable differentiation between Hepatozoon sp. 18However; a single parameter cannot be used to differentiate species using the microscopic technique.The need for an additional tool, such as the molecular technique, is needed to enable this.In addition, there is an apparent lack of relevant information regarding the genetic diversity of Hepatozoon sp. and microfilaria isolated from crocodile monitors in Thailand.The current study is the first investigation on the phylogeny of Hepatozoon and sheathed microfilaria isolated from crocodile monitors in Nakhon Pathom, Thailand.
The findings of this study show that crocodile monitors from Thailand can harbor various parasites, including both hemogregarines and filarial worms.The phylogenetic results revealed that Hepatozoon sp. in crocodile monitors had 99% similarity with H. caimani, and the Hepatozoon 18S rRNA gene was grouped in the same clade as crocodiles, reptiles, rodents, and amphibians.In this case, it is not possible for us to make a definitive determination about how the transmission occurs.However, it is plausible that transmission could occur either through prey-predator interactions or via vectors.This has been observed in African reptiles 19 and described in the case of H. domerguei infection in native reptiles from Madagascar 20 .Prey-predator transmission occurs when a predator ingests infectious cysts present in its prey.Additionally, to confirm vector-borne transmission, it is necessary to ascertain and identify the developmental stages of arthropod vectors.
In this study, the oligonucleotide pairs HepF300/900 and HEMO1/HEMO2 were used to amplify the Hepatozoon 18S rRNA gene 21,22 .This method had already been successful in inferring phylogenetic relationships between Hepatozoon spp.from snakes 23,24 .Oligonucleotides 18S and 5.8S have also been used to successfully infer phylogenetic relationships between Hepatozoon spp.from reptiles, amphibians and mammals [23][24][25][26] .The utilization of oligonucleotides HEMO1 and HEMO2 enabled the identification of a new species of Hepatozoon in Coluber constrictor priapus and Thamnophis sauritus sackenii.Through their application, the researchers were able to establish the phylogenetic relationship among hemogregarine isolates originating from Florida 27 .Therefore, the www.nature.com/scientificreports/18S rRNA sequence is useful for characterization and comparing phylogenetic relationships generic affiliations without prior knowledge of the sporogonic development of parasites 28 .
The phylogenetic results for the sheathed microfilaria COX1 gene in crocodile monitors revealed that it closely related to B. timori in humans from Japan and Indonesia.However, little is known about phylogenetic relationships of filarial worms in crocodile monitors.The phylogenetic relationship of filarial worms in wild endemic reptiles from Madagascar has been previously described 20 .In Malaysia, a new genus, Malayfilaria, along with a new species, M. sofiani, was identified in common tree shrews using the COX1 and 12S rRNA genes and the ITS1 region 29 .M. sofiani appears most closely related to Wuchereria spp.and Brugia spp.; however, it differs in several morphological characteristics.Therefore, it is important to assess the real prevalence of this parasite and investigate its implication for the host, as filarial nematodes, such as B. malayi and B. timori, are known to cause lymphatic diseases in humans living in tropical areas, while B. pahangi infects carnivores and causes zoonotic diseases in humans 29 .Unfortunately, in the current study, the PCRs targeting the 18S and the 12S rRNA genes failed to amplify sheathed microfilaria DNA.However, parasitemia levels may have been lower, thereby potentially resulting in the failure to detect parasites using PCR.Therefore, the infection levels reported in previous studies that used different primers should be compared with caution.

Conclusions
This report presents the first findings on the molecular detection of Hepatozoon sp.18S rRNA gene and sheathed microfilaria COX1 gene in crocodile monitors from Thailand.The results from the molecular analysis indicate that the evolutionary distance between the Hepatozoon sp.18S rRNA gene and sheathed microfilaria COX1 gene is greater than the distance between the previously known species, H. caimani and B. timori, respectively.Consequently, further research focusing on the transmission, interactions between hosts and parasites, and distribution of vectors for these parasites is of utmost importance, particularly in crocodile monitors.

Collection of blood samples and morphological study of the parasites
Two crocodile monitors (V.salvator) were captured in Nakhon Pathom province and restrained using a nose pole before being transported to a veterinary hospital.A veterinarian collected peripheral blood samples (n = 2) from the caudal tail vein of the crocodile monitors by using an 18-gauge needle and transferred them into EDTA-treated tubes.The samples were collected in accordance with applicable local guidelines.The blood samples were then submitted to the Vet Central Lab.The blood samples were used in blood smears for microscopic examination.Slides were air-dried, fixed with methanol, and stained with Giemsa 30 .Giemsa-stained thin blood smears were examined microscopically to assess the presence of Hepatozoon gamonts and microfilariae, as well as erythrocyte changes caused by the presence of parasites.To examine the intraerythrocytic parasite stages, digital images were obtained and measured using an Olympus CX31 biological microscope (Olympus, Japan) at 100× magnification.The measurements, in micrometers (µm), included the length and width of the parasite, with corresponding mean and standard deviation values (mean ± standard deviation).The remaining EDTA blood sample was preserved at − 20 °C for subsequent molecular analysis.

DNA extraction, amplification, and sequencing
The collected blood samples were used for DNA extraction.DNA samples were extracted using a genomic DNA blood kit (NucleoSpin ® Blood, MACHEREY-NAGEL, Germany).This process was carried out according to the manufacturer's instructions.Initially, the detection of Hepatozoon sp.DNA was achieved using primers targeting part of the 18S rRNA gene, namely HepF300, and HepR900.Positive samples then used the primers HEMO1 and HEMO2 to amplify a partially overlapping fragment of the 18S rRNA gene to obtain a longer gene portion, as shown in Table 1.Microfilariae were detected in blood smears, and three pairs of primers were used     was used as the standard for determining the molecular mass of the PCR products.The reaction products were purified using a PCR clean-up gel extraction kit (NucleoSpin ® Gel and PCR Clean-up, MACHEREYNAGEL, Germany).Purified amplified DNA fragments were submitted for sequencing using Barcode Taq (BT) sequencing and used for subsequent phylogenetic and haplotype diversity analyses.

Phylogenetic analysis
Phylogenetic reconstructions were based on the DNA sequence alignment of positive samples.Comparisons with sequences deposited in GenBank used the nucleotide BLAST.The sequences were aligned with sequences published in GenBank using the Clustal W algorithm, available in the MEGA software, version 11.0.13 31 .Phylogenetic relationships were inferred using the Maximum Likelihood (ML) and the Bayesian inference (BI) methods in MrBayes, version 3.1.2 32.The reliability of inferred phylogenetic relationships was evaluated by the statistical calculation of 1000 replicates using the bootstrapping method 33 .A Bayesian Markov Chain Monte Carlo analysis was conducted with four Markov chains (three heated chains and one cold) for 50,000,000 generations, with the trees sampled every 1000 generations.The first 50% of the trees were discarded and the remaining samples were used to construct a Bayesian consensus tree and to infer the posterior probability.Genetic distances were assessed using distance matrices under the assumption of pairwise-distance 34 and using the Kimura 2-parameter method 35 .Similarities were evaluated using the sequence identity matrix tool in the BioEdit program, version 7.0.5.3 36 .

Haplotype diversity analysis
The DNA polymorphisms and haplotype information of Hepatozoon sp. and microfilariae sequences were determined using the DnaSP software, version 5.10.01 37 .Haplotype networks were established using the TCS network tool in the Population Analysis with Reticulate Trees (PopART) software 38,39 .

Ethical approval
This research project was approved by the Biosafety Committee of Chulalongkorn University, Faculty of Veterinary Science (IBC 2231037).The authors would like to confirm that the samples were collected in accordance with applicable local guidelines. https://doi.org/10.1038/s41598-024-54276-6

Figure 1 .
Figure 1.Giemsa staining of hemogragarine gamonts and sheathed microfilaria infections from crocodile monitors in Thailand.Black arrows indicate the presence of hemogregarine gamonts of Hepatozoon (A, B).Red arrows indicate the presence of sheathed microfilaria (C).

Figure 2 .
Figure 2. The phylogenetic tree constructed using the Maximum Likelihood (ML) and Bayesian inference (BI) method based on Hepatozoon 18S rRNA sequences (917 bp).The isolates Adelina dimidiate (DQ096835) and Adelina grylli (DQ096836) were used as an out-group and the sequence of this study is in bold.

Figure 3 .
Figure 3.The phylogenetic tree constructed using the maximum likelihood method and sheathed microfilaria COX1 sequences.The bootstrap values are shown at branching points and the sequence of this study is in bold.

Figure 4 .
Figure 4. Templeton, Crandall, and Sing network of haplotypes based on the Hepatozoon 18S rRNA gene sequences examined in Thailand and globally.The small traits between haplotypes indicate the occurrence of mutations.

Figure 5 .
Figure 5. Templeton, Crandall, and Sing network of haplotypes based on the sheathed microfilaria COX1 gene sequences examined in Thailand and globally.The small traits between haplotypes indicate the occurrence of mutations.

Table 1 .
The target genes and primers (sequence and length) used to detect and characterize crocodile monitor parasites.

Table 2 .
The similarity of the Hepatozoon 18S rRNA gene sequence from a crocodile monitor in Thailand compared with global sequences.

Table 3 .
The similarity of the sheathed microfilaria COX1 gene sequence from a crocodile monitor in Thailand compared with global sequences.

Table 4 .
Polymorphisms and genetic diversity of Hepatozoon 18S rRNA and sheathed microfilaria COX1 sequences from a crocodile monitor in Thailand compared with global sequences.